Comment on "Bound states of edge dislocations: The quantum dipole problem in two 

dimensions" 
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We show that it is possible to improve some of the numerical results contained in a recent paper 
with an optimal implementation of the methods used in that paper. A careful analysis done using 
the Rayleigh-Ritz method provides a rigorous upper bound for the energy of the ground state of an 
electron in a two dimensional potential generated by the edge dislocation, as well as precise values 
for the excited states. The extrapolation of the results corresponding to different subspaces is used 
to obtain an alternative estimate of the fundamental energy of the model. The energies of the first 
500 states that we have calculated are in perfect agreement with the expected asymptotic behavior. 

PACS numbers: 67.80.B-, 02.60.-x 
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The presence of edge dislocations in a solid may have 
strong effects on the mechanical, transport, elastic and 
superconducting properties of the solid^. The knowledge 
of the spectrum of the bound states at the edge is there- 
fore helpful to assess the change in the properties of the 
material. 

With this motivation, the authors of a recent paper, 
ref.^, have considered the Schrodinger equation (SE) for 
an electron confined in a two dimensions and subject to a 
potential V{r^ 6) = pcosd/r (corresponding to a straight 
dislocation along the z-axis): 
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cos 6 ^ _ ^ 
p * = 



(1) 



Unfortunately this equation is not exactly solvable and 
its solution requires approximate methods. In effect ref.^ 
provides a detailed list of the different approaches which 
have been applied in the literature to this problem^ri; 
the numerical estimates for the energy of the ground state 
obtained in these works are reported in Table 1 of ref.— in 
units of 2mp^/h^. The last result of the table is the im- 
proved estimate obtained in that paper discretizing the 
Schrodinger equation on a uniform square grid and work- 
ing with matrices of a maximum size of 10^ x 10^. The au- 
thors quote an error of 2% for the numerical eigenvalues 
of the planar Coulomb problem, whose exact solutions 
are known^. Unfortunately, since the precision of the 
results depends on the specific problem considered, the 
accuracy of the numerical results of refi^ is not granted 
for the dipole problem and actually one may expect on 
qualitative grounds larger errors for the quantum dipole 
problem, due to the dependence on 9 in the potential. 

In both cases, the long range nature of the potential 
and its singular behavior at r = would better be taken 
into account using a nonuniform two-dimensional grid 
(see refs.— >^ for a discussion of collocation methods with 
nonuniform grids). 

A second difficulty of the computational approach 
based on the discretization of the Schrodinger equation, 
which is mentioned in refj^, is the limited number of 
states which can be obtained with acceptable precision 



with a given grid, due to the different length scales of the 
excited states: as a result the estimates for the first few 
states of the dipole potential cannot be equally accurate. 
Finally, one should bear in mind that the discretization 
of the SE, however accurate it could be, does not provide 
upper bounds on the energy of the fundamental mode: 
the result obtained in ref.^ may thus either over or un- 
derestimate the exact energy. 

The second approach discussed by the authors of ref.— 
is what they call a " Coulomb basis method" , which is es- 
sentially a Rayleigh-Ritz (RR) approach which uses the 
basis of the planar Coulomb problem (see ref.—). The 
Rayleigh-Ritz approach, in contrast to the real-space di- 
agonalization method mentioned earlier, does provide an 
upper bound to the energy of the fundamental state and 
a direct decomposition of the approximate solutions in 
the basis chosen. Unfortunately the potentialities of this 
method have not been fully exploited in that paper: the 
main purpose of this Comment is then to obtain stricter 
bounds for the energy of the ground state working with 
a larger set of functions. 

We proceed to describe the approach; the 2d hydrogen 
wave functions are^ 



with n = 1, 2, . . . and —n + l<l<n — I. 

The angular and radial parts respectively read^^ 
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1 <l <n-l 

1 = 
-n + 1 < I < -1 



(3) 



and 



{n + 
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2|Z|!y (2n- l)(n- |;| - 1)! 
• e-P"^/^ iFii-n+\l\ + l,2\l\ + l,/3nr) (4) 

where /3„ = (2k-i) ^1^- Here iFi{a,b,c) is the confluent 
hypergeometric function. 
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The bound state energies of the 2D hydrogen are sim- 
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It is convenient to introduce the notation a = and 
regard a as a variational parameter controlhng the length 
scale. We then rewrite the matrix elements of the Hamil- 
tonian making the dependence upon a explicit: 



The RR approach requires the calculation of the matrix 
elements of the Hamiltonian of eq. ([1]) in the basis ([2]): 



{nili\H\n2l2) = Si^ijn^ ^ri2^ni 



RniiAr)Rn^i^{r)dr 



where 



r(2) 



Xh{0) cosexb_iO)d0 = [Si^+i^i., +Si^^i^i^] 
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The evaluation of the radial integral is not straightfor- 
ward as in the case of the angular integrals but it can 
also be done analytically. Defining = —ni + \li\ + l 
and bi = 2\li\ + 1 { i = 1,2), the confluent hypergeo- 
metric functions in this integral reduce to polynomials of 
degrees —ai and —02 respectively, for ai,2 < and 61,2 
positive integers. The original integral is therefore re- 
duced to a sum of integrals which can be done explicitly: 



Rmh{r)Rn2i2{r)dr 
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where the {a)k are the Pochhammer symbols and 

^nuh,n;n.M,n = 4(2n2-l)l'il+^-i-3(2ni-l)l'^l+^-^- 
(-ni + |;i| + l),,(-n2-f |;2| + l)j. 
• r(ji+j2 + |;i| + |^2| + l) 
^jT{ni + \h\)T{n2 + \h\) 

= r(ji + i)r02 + 1) 
r(ji + 2|Zi| + i)r(j2 + 2|Z2| 



(7) 
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The authors of ref.^ have used the analytical expres- 
sions for these integrals, although they had to restrict 
their calculation to a set of 400 basis functions, corre- 
sponding to — n + 1 < Z < n — 1 and 1 < n < 20, due 
to the numerical round-off errors that become important 
for the matrix elements corresponding to larger quan- 
tum numbers. In our numerical calculation we have used 
Mathematica 8—, obtaining symbolic expressions for the 
matrix elements, which were then evaluated numerically 
avoiding the round-off errors which would appear in a 
fully numerical calculation. 
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This is precisely the approach followed by Dasbiswas 
et al. in ref.^, who observe that the bound for the ground 
state energy obtained for a = 1 is not good: after mini- 
mizing with respect to a they obtain an improved bound 
El > —0.1257 working with 400 states (we adopt their 
convention of reporting the energies in units of 2mp'^ /h^). 

Interestingly, these authors use the original basis {a ~ 
1) for the excited states, claiming that "the real-space 
diagonalization methods provide a better estimate of the 
ground-state energy whereas the Coulomb basis method 
is more suitable for higher excited states. " It is not clear 
on what grounds this statement is made and actually we 
will show in this Comment that the choice of calculating 
the excited states for a = 1 is far from optimal. We will 
also obtain an alternative estimate for the energy of the 
ground state of the quantum dipole problem, which falls 
above the one calculated in ref.^. 

Our first observation concerns the number of bound 
states which are obtained in the calculation at a given 
a: while for a = 1 there are 149 bound states, for 
c^var = 89.57 (corresponding to the minimum of Ei) 
there are just 3 bounds states. This behavior should not 
surprise us, since a determines the radial length scale 
of the wave function, and different states have different 
range: the fact that for a large a fewer bound states are 
present, simply tells us that the physical states which 
are not captured by the calculation have a larger range. 
Therefore, if one wants to estimate a few excited states, 
one needs to calculate these states using different values 
of a to account for the different length scales of each 
state. 

The fundamental question is therefore how to choose 
a: if the problem under consideration has certain sym- 
metries, which are also symmetries of the basis, then the 
variational principle applied to a trial wave function with 
that symmetry will provide again an upper bound to the 
lowest mode in that symmetry class: for our problem 
this is the case of functions which are odd with respect 
to the change y — ?> —y. For the remaining states, the 
eigenvalues of the RR matrix will vary with a, without 
providing a variational bound. However, if we consider a 
given state, its exact energy and wave function will be in- 
dependent of a, which is an unphysical parameter of the 
basis. As such, we may argue that the optimal value of a 
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FIG. 1: (Color online) Variational bound for E\ as function of 
the maximum quantum number used in the RR approach (a) 
and the extrapolation of E\ using the Shanks transformation; 
the horizontal line corresponds to the fit of eq. pup . (b). 



(and correspondingly the most accurate value for the en- 
ergy) will be the one for which the calculated eigenvalue 
is less sensitive to changes of a. This is the essence of the 
" principle of minimal sensitivity" (PMS)^^ . We will soon 
use the PMS to calculate the excited states of eq. ([T]) . 

We now proceed to illustrate our numerical results, 
concentrating for the moment on the ground state. 

In the first plot of Fig [T] we show the variational bound 
for the ground state energy of eq. (P) as a function of the 
maximum quantum number used in the RR method (we 
have considered even values of rimax)- For a given rimax, 
the portion of Hilbert space used in the calculation con- 
tains n"^^^ states (the results obtained by Dasbiswas et 
al., for instance, correspond to Umax = 20 and therefore 
to a 400 X 400 matrix). For each value of Umax we have 
selected the value of a for which the corresponding Ei is 
minimal; the numerical calculation is stopped when the 
convergence on the first 20 digits of Ei was reached. 

Remarkably these data display a very regular behavior, 
which is well described by the fit 
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-0.131678 



0.0184 



(0.661n+ 1) 



0.426 



(10) 



In the second plot of Fig IT] we show the values for 
the ground state energy obtained repeatedly using the 
Shanks transformation on the sets of the energies of the 
first plot. The index k in the horizontal axis indicates to 



the number of consecutive Shanks transformations used, 
while the value on the vertical axis is the correspond- 
ing result obtained with the Shanks transformations, in- 
volving the values corresponding to largest n„iax- These 
values seem to converge to the value Ei « —0.1314. 
The horizontal line corresponds to the constant value 
El = -0.131678 reached by eq. ^ for n^ax oo. 
Notice the disagreement of our results with the value ob- 
tained in refji, Ei ~ —0.139. Being fair, there is no rig- 
orous criterion granting that our result is more precise, 
although the extrapolation of very regular sequences of 
numbers typically provides very accurate estimates. 

The behavior of the optimal a obtained with the PMS 
at different nmax is described very well by the a cubic fit 
Q,(FiT) ^ _0.00003847n3^„^ + 0.2104n2^„^ + 0.3062?7™a:,- 
0.4347. There is a clear physical justification of the be- 
havior of a, which grows with Umax- as the number of 
states in the calculation is increased, one can use a ba- 
sis with shorter length scale (i.e. larger a) to build the 
approximate eigenfunctions of the problem. 

Although we believe that the values of Ei that we have 
obtained with the fit pUj) or with the Shanks transfor- 
mation are more precise than the result obtained in refii, 
we are aware that the extrapolation of the values of Ei 
do not themselves fulfill a variational bound. In other 
words, we may expect them to be closer to the exact 
value, but we cannot be sure that they fall above it. 
We will now show that it is possible to obtain stricter 
bounds on Ei, even working with less states. As we have 
mentioned before, the largest set of states used in FigH] 
corresponds to Umax = 60, i.e. 3600 states (1 < n < 60 
and —n -I- 1 < Z < n — 1). We may wonder to what ex- 
tent the results would change by restricting the states to 
< I < Imax- While one can safely drop the negative 
values of Z, since the ground state must be symmetric 
with respect to y — ^ —y, it is not a priori clear the error 
introduced by using an upper cutoff Imax- 

Let V = (vi,V2, ■ ■ ■ ,vzmo) be the eigenvector corre- 
sponding to the smallest eigenvalue of the 3600 x 3600 
matrix obtained in RR approach: the normalization of V 
implies that X]^=T ^"^ ~ ^- deviation from 1 of the 
sum when the components of V corresponding to I > Imax 
are set to zero will give us an idea of how important these 
states are for the calculation. For Imax — 5 one finds out 
that this deviation is completely negligible, (5 = 2.2 x 
10~^. To confirm this finding we may compare the lowest 
eigenvalue of the full 3600 x 3600 matrix, Ei < ""^ = 
-0.12788587515854182342, with the lowest eigenvalue of 
the reduced 345 x 345 matrix, corresponding to Imax — 5, 
El < e[""^'''""^'^ = -0.12788587377192016287. We have 



^(reduced) _ ^jiuu; ^ ^ 4 ^ iQ-'J, which is of the Same 
order of the deviation discussed above. The effective di- 
mensionless coupling constant calculated with the wave 
function corresponding to [g g — 0.017, which 

agrees with the result found in refii, using a simple vari- 
ational ansatz. 

With this result in mind we have built a 585 x 585 
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matrix, corresponding to Umax = 100 and Imax = 5, ob- 
taining the bound Ei < -0.12864468596110909173. No- 
tice that the fit of eq. (fTO|) for rimax = 100 provides a 
result which is very close to this, £;f' ^ = -0.128612. 
We have then built a 885 x 885 matrix, corresponding to 
nmax — 150 and Imax — 5, obtaining our most precise 
bound El < -0.1291697936750557573. Also in this case 
the fit of eq. PH)) for n„iax = 150 provides a result which 
is very close to this, e[^^^ = -0.129093. 

We now discuss the excited states of eq. ([T]): as we 
have mentioned earlier there is no valid reason for cal- 
culating the excited states at a = 1. Fig. 4 should con- 
vince the reader of this point: here the energy of the 
fifth state is plotted at different values of a, and com- 
pared with the value at a = 1 (dashed line), which 
corresponds to the choice done in^. The dotted line at 
the bottom corresponds to the minimum of the curve, 
and it provides the most accurate value for Eq which 
can be obtained working with a subspace correspond- 
ing to Umax = 60. Table U illustrates the different re- 
sults obtained using the two approaches for the first 
five states. The last column reports the quantity A„ = 
[Ek''^'^ - E^r'^)l{Ei^'''^ + E^^='\ which provides 
an estimate of the error done using a = 1. Interest- 
ingly the values obtained with the PMS for the third and 
fifth states are close to the ones obtained in refji using 
a discretization of the Schrodinger equation: these states 
correspond to smaller values of apMSi indicating that 
their length scales are larger than those of the other three 
states. So, for example, the third state, has a smaller a 
than the fourth state, which is higher in energy. The au- 
thors o£i have also observed a similar behavior for some 
of the states (23*'* and 24*'') that they have calculated, 
although they have not given "any satisfactory explana- 
tion for these irregular features" . Our explanation of this 
phenomenon is simple: if a state with a modest spatial 
extent has a larger probability density in the region close 
to a; = (recall that the potential is attractive for a; < 0), 
its energy may be higher than the energy of a state with 
larger spatial extent but smaller probability density in 
the region close to a; = 0. 

In Fig. [3] we plot the optimal value of a as a function 
of the level number n, obtained using the subspace cor- 
responding to Umax = 60. The assumption used in refi^, 
a = 1, is clearly valid only for the higher excited states 
(n > 100). Notice the oscillations of apMS, which signal 
the presence of contiguous states with different length 
scales, as mentioned earlier. 

In Fig. [4] we have compared the behavior of —1/En 
obtained using either the PMS (solid line) or setting a = 
1 (dashed line), with the asymptotic law of eq.(lO) of ref>i 
(dotted line). The plot clearly proves the superiority of 
the PMS results. This superiority can also be established 
by looking at the fit of the first 500 values oi ~1/En: 

-l/-^«L=i = 16.0528n-|- 1.384v^- 2.085, 
-l/^^nlap„s ^ 15.9897n + 0.431 V^- 9.162 . 



TABLE I: Comparison of the energies of the first five states of 
the "dipole" potential calculated with a = 1 (second column) 
and using the PMS (third column) for Umax ~ 60. Here 







— 1^ N //-,-,( P MS) 
0/(£^n ' 


r-,(a — l) \ 




n 




p(PMS) 


CtPMS 


A„ 


1 


-0.0970117 


-0.127886 


767.132 


0.137281 


2 


-0.0328379 


-0.0394579 


294.189 


0.0915679 


3 


-0.0220914 


-0.0232932 


137.674 


0.0264818 


4 


-0.016764 


-0.0193729 


161.317 


0.072193 


5 


-0.0119611 


-0.0125862 


86.4652 


0.0254668 
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FIG. 2: (color online) Approximate energy of the fifth state 
of the "dipole" potential, as a function of the variational pa- 
rameter a, for Umax = 60. The dotted line corresponds to 
the minimum value, while the dashed line corresponds to the 
value for a = 1, as in ref.— . 



In conclusion, we believe that we have presented nu- 
merical results which complement and extend the re- 
sults of reffi; we have shown that the application of the 
Rayleigh-Ritz method to this problem provides very ac- 
curate results, both exploiting the Shanks transformation 
and performing the calculation on a selected portion of 
the Hilbert space. In this way we have calculated the 
ground state energy corresponding to rimax = 150, work- 
ing with a matrix of size 885 x 885, instead of the full 
matrix 22500 x 22500, obtaining our most precise bound, 
El < —0.1291697. Using the numerical results for the 
excited states, calculated using the PMS, we have re- 
produced with high accuracy the expected asymptotic 
behavior of eq.(lO) of ref.—. 
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FIG. 4: (color online) —1/En for the excited states between 
n — 400 and n = 500. The solid line is the result obtained 
with the PMS, while the dashed line is the result obtained 
using Q = 1. The dotted line is the asymptotic law of eq. 
(10) of ref.-. 
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